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Abstract 

The control of the flutter instability and the conversion of the dangerous character of the 
flutter instability boundary into the undangerous one of a cross-sectional wing in a 
supersonic/hypersonic flow field is presented. The objective of this paper is twofold: i) to 
analyze the implications of nonlinear unsteady aerodynamics and physical nonlinearities on the 
character of the instability boundary in the presence of a control capability, and ii) to outline the 
effects played in the same respect by some important parameters of the aeroelastic system. As a 
by-product of this analysis, the implications of the active control on the linearized flutter 
behavior of the system are captured and emphasized. The bifurcation behavior of the open/closed 
loop aeroelastic system in the vicinity of the flutter boundary is studied via the use of a new 
methodology based on the Liapunov First Quantity. The expected outcome of this study is: a) to 
greatly enhance the scope and reliability of the aeroelastic analysis and design criteria of 
advanced supersonic/hypersonic flight vehicles and, b) provide a theoretical basis for the 
analysis of more complex nonlinear aeroelastic systems. 
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Introduction 


During the evolution of the combat aircraft, dramatic drops of the flutter speed can be 
experienced. As a result, the aircraft can attain the critical flutter speed and depending on the 
nature of the flutter boundary, i.e. catastrophic or benign, the aircraft can exhibit dramatic 
failures or can survive, respectively. At this stage, according to the flight regulations, the flutter 
speed should be (22 - 25)% larger than the maximum speed the airplane can experience in the 
dive flight. This large margin of security is imposed, in order to prevent the catastrophic failure 
of the aircraft, in the event that the flutter speed would be crossed, and based on what was 
considered until now that its crossing would result in a catastrophic failure of the aircraft. 
However, as it was shown, the flutter boundary can also be benign, when, in this case, it can be 
crossed without the occurrence of catastrophic failures. 

This suggests the considerable importance of determining the conditions that result in the benign 
and catastrophic characters of the flutter boundary, and in the development of proper 
mechanisms able to convert, in an automatic way, the catastrophic flutter boundary into a benign 
one. These facts emphasize the considerable importance of at least two issues: a) of including in 
the aeroelastic analysis the various nonlinear effects, on which basis is possible to get a perfect 
understanding of their implications upon the character of the flutter boundary, and b) the 
importance of implementing adequate active feedback control methodologies enabling one not 
only to increase the flutter speed, but convert the catastrophic flutter boundary into a benign one. 
The problem of aeroelastic stability in the vicinity of the flutter boundary, for both the lifting 
surfaces and supersonic panels, is an issue that requires a great deal of research toward 
understanding the effects of the nonlinear unsteady aerodynamics and structural/physical 
nonlinearities. Related to this problem, as it was shown in Ref. 1, the various nonlinearities that 
occur in the aeroelastic governing equations and are of structural (i.e. arising from the 
kinematical equations 2 ' 7 ), physical (arising from the constitutive equations), or coming from the 
unsteady aerodynamic loads 2-8 9 , can render the flutter boundary either catastrophic or benign. In 
the former case, crossing the flutter boundary results in an explosive failure of the structure, 
while in the latter case, a monotonous increase of displacements with the increase of the flight 
speed takes place. 
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It is well known that, aerodynamic and structural nonlinearities affect the aeroelastic response of 
the wing and the flutter boundary 1 ' 3,8 ' 9 . Taking in to account these nonlinearities, an 
understanding of their potential influence on the character of the flutter can be reached. 

This paper restricted to a supersonic/hypersonic cross-sectional wing is intended to address these 
issues. 

Due to the practical importance of these problems, a thorough investigation is needed not only to 
determine whether or not, for the actual geometrical physical and aerodynamical parameters of 
the structure, the flutter boundary is of a catastrophic nature, but also to be able to convert it, if it 
is the case, to a benign one. 

In the case when, due to the character of considered non linearities coupled with that of geometric 
and physical parameters of the lifting surface, the flutter boundary is catastrophic, an active 
control mechanism should be implemented as to convert the flutter boundary into a benign one. 
In this sense, by controlling the aeroelastic response, the active control methodology 31015 is 
likely to play a powerful role toward avoiding the occurrence of catastrophic failures, toward 
rendering the flutter boundary a benign one, and the expansion of the flight envelope. 

In this paper a general method able to approach both the problem of lifting surfaces flying at 
supersonic/hypersonic flight speed regimes and that one of panel flutter 1 , have been developed. 
This methodology will permit to identify the nature of the flutter boundary (i.e. benign or 
catastrophic) that, in contrast to the actual brute methods, is of an analytical nature and include 
all the factors, starting with the nonlinearities and continuing with all the aeroelastic parameters 
of the structure. In this context, general conclusions related to this matter can be reached. 
Moreover, the analytical methodology to address this issue is based on powerful mathematical 
tools elaborated by famous mathematicians, such as Bautin 16 , a disciple of Liapunov, Hopf, etc. 
All these elements suggest that, in order to predict the flutter instability boundary and its 
character, it is vital to efficiently model the unsteady aerodynamics of lifting surfaces associated 
with the corresponding flight conditions. 

It should be noticed, that in spite of the great practical importance, the literature dealing with the 
problem of the determination of the flutter boundary of sectional-wing and its character in the 
presence of both physical and aerodynamic nonlinearities in the conditions of the 
supersonic/hypersonic flight speed in the subcritical flight conditions is quite void of any such 
results 1 . 
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There are many potential sources of nonlinearities, which can have a significant effect on an 
aircraft’s aeroelastic response. One essential limitation involving the linearized analysis is that it 
can only provide information restricted to the flight speed at which the aeroelastic instability 
occurs. Furthermore, the linearized analyses are restricted to cases where the aeroelastic response 
amplitudes are small. Often this assumption is violated prior to the onset of instability. Thus, to 
study the behavior of aeroelastic systems in either the post-instability region or near the point of 
instability, the nonlinearities of a physical geometrical or aerodynamic nature must be accounted 
for. 

Hopf-Bifurcation 

The issue of the character of the flutter boundary, i.e. benign or catastrophic , can be revealed 
via determination of the nature of the Hopf-Bifurcation (i.e. supercritical or subcritical, 
respectively) as featured by the nonlinear aeroelastic system 1,17 . 

A pictorial representation of the behavior of the aeroelastic system in the vicinity of the flutter 
boundary versus the dimensionless flight speed in terms of the growth of the limit cycles 
oscillation (LCO) amplitudes is depicted in Figs, la and lb. The curves indicated by 1 and 2 are 
characteristic to the stable and unstable domains. In these two cases, below the linear flutter 
speed V F , the aeroelastic system is stable, whereas above V F , the system is unstable for any 
value of the initial conditions. Above V F in 3 the system tends to the stable LCOs, that is known 
as supercritical Hopf-bifurcation (H-B), identifying an undangerous flutter boundary, whereas 
the curve indicated by 4 characterizes a subcritical Hopf- Bifurcation that yields the catastrophic 
failure of the spacecraft structure. The latter character of the flutter boundary is referred to as 
dangerous. The solid lines 3 and 4 are also known as pitchfork bifurcations, whereas the dotted 
lines identifies the two cases known as knee-like shape bifurcations in which a turning point is 
experienced. Moreover, according to the initial conditions, in 3 and 4 the LCOs amplitudes of the 
system decrease or increase (stable, unstable). In addition, the cases depicted in 5 and 6 are 
stable and unstable. Due to its highly destructive effects, flutter instability must not occur within 
and on the boundary of the operating envelope of the flight vehicle. In this sense, the nonlinear 
approach of lifting surfaces of aeronautical and space vehicles permits determination of the 
conditions under which undamped oscillations can occur at velocities below the flutter speed, (in 
this case the flutter boundary is catastrophic) (case corresponding to curve 3 in Fig. 1 ), and also 
of the conditions under which the flight speed can be exceeded beyond the flutter instability, 
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without catastrophic failure (in this case the flutter boundary is benign ) (curve 4 in Fig. 1). As a 
by-product of the nonlinear aeroelastic analysis, the flutter speed V F i.e. the speed for which the 
undisturbed form of the considered structure ceases to be stable, can be determined via a linear 
stability analysis of the aeroelastic system. 

Loosely speaking, the Hopf-bifurcation theorem stipulates that if the characteristic equations of 
the linearized system about an equilibrium position exhibit pairs of complex conjugate 
eigenvalues that cross the imaginary axis as one of the control parameter varies, (in the present 
case this parameter is the flight speed V ), then for the near-critical values of V there are limit 
cycles close to the equilibrium point. Just how near the critical flutter speed V F has to be is not 
determined, and unless a certain rather complicated expression is nonzero, the existence of these 
limit cycles is only assured exactly at the flutter speed (y =V F ). The sign of the expression 
determines the stability of the limit cycle, and whether the limit cycle exists for subcritical 
(F > V r ) or supercritical (V < V r ) parameter values. 

Behavior in the Vicinity of the 
Flutter Instability Boundary 

In this study, the determination of the dangerous/undangerous character of the critical flutter 
boundary and its control reduces to the determination of the sign of the first Liapunov quantity, 
for the instability boundaries that correspond to the purely imaginary roots of the characteristic 
equation. The approach used in Refs. 1 and 2 for the panel flutter, is extended in the present 
work to the aeroelasticity of nonlinear 2-D lifting surfaces. 

The behavior of the general dynamic systems near the boundaries of the stability domains 
was investigated by Bautin 16 who considered only those portion of boundary of the region of 
stability for which the characteristic equations exhibits either one root only, or two roots, that are 
purely imaginary. 

Starting with the expression of the first Liapunov quantity L(V f ), pertinent conditions defining 
the dangerous and undangerous character of the flutter critical boundary can be determined on 
the basis of the considered nonlinear aeroelastic system. Specifically, from the condition that 
L(V f ) should satisfy the inequalities L(V f )<0 and L(V / y; )>0, we determine the undangerous 
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and dangerous portions, respectively, of the boundary of the Routh-Hurwitz domain, and 
implicitly, the influence played by the various nonlinearities included in the system. 

In the next developments the first Liapunov quantity L(V f ) corresponding to the nonlinear 
flutter equations of open/closed loop cross-sectional wing in a supersonic/hypersonic flow field 
is derived and is used toward determining the conditions that characterize the nature of the flutter 
stability boundary. 

Nonlinear Model of Cross-Sectional Wing 
Incorporating an Active Control Capability 

Toward formulating the open/closed loop aeroelastic theory of cross-sectional wing in a 
hypersonic flow field both the aerodynamic and physical nonlinearities are included. This is 
motivated by the fact that these nonlinearities can contribute differently to the character of the 
flutter boundary. Moreover in the case when, due to the character of considered nonlinearities, 
the flutter boundary is catastrophic, an active control mechanism will be implemented as to 
convert the flutter boundary into a benign one. 

The open/closed loop aeroelastic governing equations of a cross-sectional wing featuring 
plunging and twisting degrees of freedom, elastically constrained by a linear translational spring 
and nonlinear torsional spring, exposed to a supersonic/hypersonic flow field are 19 ' 20 : 

mh(t)+ S a a(t)+c h h(t)+ K h h(t)= L a (t), (1) 

S a h(t)+ I a d(t)+ c a a(t)+M a =M a (t)-M c . (2) 

Herein h(t) is the plunging displacement (positive downward), a(t ) is the pitch angle (positive 
nose up), the superposed dots denote differentiation with respect to time t, m is the structural 
mass per unit span, S a is the static unbalance about the elastic axis, I a is the mass moment of 

inertia about the elastic axis of the airfoil, c h and c a , K h are the linear plunging and pitching 
damping and stiffness coefficients, respectively. Moreover, in Eq. 2: 

M a =K a a{t)±8 s K a ^{t), (3) 

represents the overall nonlinear restoring moment, that includes both the linear and the nonlinear 
stiffness coefficients, K a and K a , respectively. 
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Within a linear model ( 8 S =0), M a would simply be replaced by K a a(t). The nonlinear 
coefficient in Eq. (3) can assume positive or negative values. In the former case we deal with a 
hard physical nonlinearities, while in the latter one with a soft physical ones. Notice that this 
nonlinearity appears in the present case in the equation relating the restoring moment with the 
pitch angle and, it has as such the character of a constitutive equation. For this reason, in contrast 
to the terminology of structural nonlinearity, that in the authors’ opinion is not proper for this 
case, herein we are using that of physical nonlinearity. 

The tracer 8 S that identifies this type of nonlinearities can take the values 1 or 0 depending on 
whether this nonlinearity is accounted for or discarded, respectively. 

The active nonlinear control can be represented in terms of the moment M c as: 

M c = M c a(t)+8 c M c a*(t), (4) 

where M C ,M C are the linear and nonlinear control gains, respectively. Also in this case, within 
a linear active control methodology the tracer assumes the value 8 C = 0 . 

As to reduce the aeroelastic governing equations to a dimensionless form, we also define the 
parameter: 

B = kJk* , (5) 

and the two active control gain parameters: 

Vi=M c /K a , (6a) 

=M C /K a . (6b) 

B constitutes a measure of the degree of the nonlinearity of the system and corresponding to 
B< 0 or £>0, the nonlinearities are soft or hard, respectively, while denote the 

normalized linear and nonlinear gains, respectively. 

Let now derive the nonlinear unsteady aerodynamic lift and moment. From Piston Theory 
Aerodynamics (PTA) 21 ' 2 ' 3 , the pressure on the upper and lower faces of the lifting surface can be 
expressed as: 

( . y>7(y-i) 

p(x,t) = 1 + ^— — , ( 7 ) 

2 a „ 

\ ~ / 

where the downwash velocity normal to the lifting surface v z is expressed as 1 : 
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and 


v, = 



0W 

0.V 


\ 

sgnz , 

) 



(8) 


(9) 


Herein p^,p x , t/„ and a x are the pressure, the air density, the air speed, and the speed of sound 
of the undisturbed flow, respectively. In addition, K is the isentropic gas coefficient, ratio of 
specific heats, under constant pressure and constant volume, respectively ( /c = 1 .4 for air), w is 
the transversal displacement of the elastic surface from its undisturbed state 
w(t) = h(t) + a{tXx- x EA ) y and sgn z , assumes the values 1 or -1 for z > 0 and z < 0 , respectively. 
In addition, x EA ~bx 0 is the streamwise position of the pitch axis measured from the leading 
edge (positive aft), whereas b is the half-chord length of the airfoil. 


Retaining, in the binomial expansions of (7), the terms up to and including (v 7 /a„,) 3 , yields the 
pressure formula for the PTA in the third-order approximation 2 ' 22 : 


P , v - 
— = 1 + K—y + 


tc(k: + 1) 
4 



! «fr + !) 
12 



( 10 ) 


Herein the aerodynamic correction factor y 23 : 



enables one to extend the validity of the PTA to the entire low supersonic-hypersonic flight 
speed regime. As mentioned in Ref. 1, Eq. (10) are satisfactory even for M >2. It is also 
important to remark that, Eqs. (7) through (10) are applicable as long as the transformations 
through compression and expansion may be considered as isentropic, i.e. as long as the shock 
losses would be insignificant (low intensity waves). 

On the other hand, a more general formula for the pressure, obtained from the theory of oblique 
shock waves (SWT) 24 , and valid over the entire supersonic - hypersonic range, can be applied 
(see Ref. 2, 22, 24, 25): 


_P_ 

P~ 


, v. k(k + l) 

[ + k— y + — 

a 4 



, gfr + 0 2 

32 



( 12 ) 
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Eqs. (10) and (12) are similar, the only difference occurring in the cubic term. This is explained 
by the fact that the entropy variation appears in the pressure expansion, beginning from the third- 
degree terms. In contrast to the previously mentioned formulae, this one, holds valid both for 
compression prior the shock wave and for expansion (it is being assumed that the shock waves 
are all the time attached to the sharp leading edge, and that the flow behind these waves remains 
supersonic). Moreover, Eq. (12) encompasses a number of advantages such as: i) takes into 
account for the shock losses occurring in the case of strong waves, ii) can be used, in the form 
derived in Refs. 24, 25, over a larger range of angles of attack (a <30°) and Mach numbers 
(M >1.3), and iii) is still valid for Newtonian speed (M — »°°;y — »l). Comparisons of results 
depicting the dangerous and undangerous character of the flutter boundary using these two 
formulae will be presented at the end of this paper. However, as it clearly appears, within the 
linear stability analysis, the flutter speed evaluated with these two expressions does not exhibit 
any differences. This is due to the fact that for such an analysis, only the linear terms are retained 
toward evaluation of the flutter speed, and these ones are identical for both PTA and SWT (see 
Eqs. (10) and (12)). 

The two cubic term coefficients differs by 10% for k = 1.4, and so, for a better prediction of the 
character of the flutter instability boundary, it has to be included (see Ref. 22). Moreover, the 
SWT, which considers the pressure losses through shock waves, gives more accurate results then 
the PTA. 

The flutter speed vs. flight Mach number obtained from the PTA and SWT including (i.e. 
extending the flight speed toward low supersonic range), and discarding the correction factor for 
the speed y (in the sense of considering M l » 1 , and consequently y — » 1) is shown in Fig. 2. 
In addition, the flutter boundary obtained via the use of the exact supersonic unsteady 
aerodynamics is depicted on the same figure. In the low supersonic flight speed regime, the PTA 
and SWT with the corrective term gives a rather good agreement with the supersonic flow theory 
and so, this correction should be included. At the same time, for higher supersonic Mach flight 
numbers, the differences in the flutter predictions based upon the indicated aerodynamic theories 
practically disappear. In the next developments, unless otherwise stated, Piston theory 
aerodynamics has been applied. 
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Considering that the flow takes place on both surfaces of the airfoil (with the speed 
U* = U~ = U „ ), from the Eqs. (7) - ( 10) we can express the aerodynamic pressure difference dp 
as: 


8p\ 


PTA 





+ VV, , 


K 


l + K 2 ,,■>/ 1 

7 M^(w, — +w, 

12 ' U„ 



(13) 


Herein, M „ =U„/a x is the undisturbed flight Mach number, while q = -j pJJt is the dynamic 


pressure. 

Notice that, in Eq. (13) and in the next developments, the cubic nonlinear aerodynamic term will 
be reduced to the (w, t ) only, in the sense that the nonlinear aerodynamic damping will be 

discarded. There is no clear-cut for the start of the hypersonic flight speed regime. Generally 
speaking, speeds above Mach-5 are considered hypersonic. This is the speed at which the 
aerodynamic heating becomes important in aircraft design. Having in view that this effect will 
not affect the conclusions about the implications of the considered nonlinearities, the effects of 
the non-uniform temperature field are neglected here. 

Next, the nonlinear unsteady aerodynamic lift L 0 (/) and moment M a {t) per unit wing span can 

be obtained from the integration of the pressure difference on the upper and lower surfaces of the 
airfoil: 

Ai(0= and M a {t) = -j 2 o b 8p{x-x EA )dx. (14a,b) 

Their nonlinear final expressions can be cast in compact form as: 

L. Mlujl + r)yW(t )+l2[>,(l)+ (b - x a M<)]}, (15a) 

M .(')=^khr{\ 2 v.(b-x E M)+Miu.(b-x CA t\+Ky-a > (,) 



+ 4(3(6 - )h(t)+ (4 b 2 - 6 bXu + 3.4, )*(0]} ■ ( 1 5b) 

As a result, the governing equations (1) - (2) considered in conjunction with Eqs. (15a,b) feature 
the inertial and the aerodynamical types of couplings. 

Upon denoting the dimensionless time T = Uj/b, the system of governing equations, including 
the control terms, can be cast as: 

£ '(t)+ *„<*'(*)+ 2C„ {a/v )§'(t)+ (m/v ) 2 £(t)= l a (r), ( 1 6) 
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(Xa hi )r(*)+ <*'(*)+ ( 2 C« /V )*'(*)+ I/V 2 a(r)+j/v 2 _Ba 3 _(j_)= 

= m i) (T)- y l /V 2 «(T)-y 2 /V 2 « 2 (T) > (17) 

in which l u (= L a b/mUl ) and m a (= M a b 2 / 1 a U l), denote the dimensionless aerodynamic lift and 
moment, respectively. These are expressed as: 

L (?) = - ^ {l2«(r)+ Mi (1 + tc)y V (t)+12[£'(t)+ (ft - )/fta'(r)]}, ( 1 8a) 

12^M m 

m “ ^ = vT^aa — t! 1 2 ( b ~ x EA )/ ba ( T )+ M l {b - x EA )/b (1 + x)y 2 a 2 (r) 

127 M/ m r a 

+ 4[(3(ft-^)/ftr(T)+(4ft 2 -6 bx EA + 3*^ )/ft 2 «'(*))]} • 08b) 

Herein, £ = ft/ft, is the dimensionless plunging displacement, the superposed primes denote 
differentiation with respect to time T, ji-m/Apb 2 is the dimensionless mass parameter, 
= S a /mb is the dimensionless static unbalance about the elastic axis E, r a = Jl a / 2 is the 
dimensionless radius of gyration with respect to E, V =U, x ,/b(O a is the non-dimensional airspeed 
parameter, (O-(D h /03 a is the plunging-pitching frequency ratio, where 0, (=c,j2m(o h ), 
£« (= c a/ 2/ a <y a ) are the damping ratios in plunging and pitching, and where co a = JKjm and 

c° a = JkjT« are the uncoupled frequencies of the linearized aeroelastic system counterpart in 

plunging and pitching, respectively. In addition, in the above equations the terms underscored by 
a single solid and a dashed line identify the aerodynamic and physical nonlinearities, whereas the 
terms underscored by a double solid line identify the active control terms, respectively. 

As previously stated, the conditions of dangerous or undangerous character of the flutter 
instability boundary were obtained by the use of the Liapunov first quantity L(V f ). This quantity 
will be evaluated next. 

To this end, the system of governing equations is converted to a system of four differential 
equations in the form 1 ' 216 : 

{Jx 4 

— L =^a ( J ) x m +P l (x l ,x 2 ,x„x 4 ), 7 = 1,4. (19) 

m=l 

For the present case the functions P J (x l ,x 2 ,x^,x 4 ) include both the physical, aerodynamic and 
nonlinear control terms that can be cast as: 
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P ; (x p x 2 ,^,xJ=X4 y V +2 S4 ;) ^^ + £4/ ) *.- +3 ^ajpxfx, +6 ^a ( £x,x,x k .(20) 

1=1 (./=! i=l /,/=] ij,k=l 

(i*l) ' ' 


i.l= 

(i*D 


i.l,k=l 
(i*l*k ) 


Specializing these expressions for our case, the governing equations in state space form can be 
reduced to: 

dx x 

d T 
c/v. 


■ = *i. 


</r 


= x, , 


(21a) 

(21b) 


— a?^x, (t)+ a^x 2 (t)+ S A a2,\xl (r)+ <5 5 a£^x 2 (r)+ 5 c a2i 2 x\ (t )+ af^x, ( t)+ a 4 ^x 4 (t) , (21c) 

dr 


<^x 4 
d T 


= fl^X, (t)+ a ^ X 2 (t)+ ^4^222 *2 ( T )+ ^s a 222 X 2 ( 0 + ^c a 222 X 2 ( T )+ X ] (0 + O ^ X a (t) . (2 1 d) 


where (£ = x , \a = x 2 ;£' = x 3 ,cl = x 4 ). The linear active control is included in the coefficients 
af^ and a^\ whereas the nonlinear ones are included in the terms accompanying the tracer 8 C . 

The coefficients of Eqs. (21) are displayed in Appendix A. As it clearly appears from these 
coefficients, if the streamwise position of the pitch axis coincides with that of the center chord, 
the nonlinear aerodynamic terms decay in importance, in the sense that the aerodynamic pitching 
moment vanishes. 

Considering the solution of Eqs. (21) under the form x j = A j exp(fttf), the characteristic equation 
corresponding to the linearized system counterpart is: 

ft) 4 + poo* 4- q(Q 2 + r(0 + 5 = 0, (22) 


where: 

vA 4 + 34 + 6*0 (*„ - 1)- 6X„ 1+ 3r„* (Vy + 2 + 2 ) 

P ~ 3M.HVti-xl) ’ 

= + 0) 2 +y/,)+2C g (Vy + 2M 0 . j U6)Cj] 

q 3 MltfV'Xrl-zl) 

Vy[v (y + 3 Mji - 3 M„nx a )+ SMjwiC, + 6Mjmx 0 2 C, - 3 M^ 0 (V + 45^, )] > 

3 M 2 7rV 2 (r a 2 -* 2 ) 

_ 5T 2 [Mji{ 1 + yr, )r g 2 + V 2 y(x 0 - l)j 

M m fiV 4 (rZ-x 2 a ) ’ 


(23a) 


(23b) 

(23c) 
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= 3r„ 2 [2M^(W 2 C a +{\ + y/ f XVy + 2M^G)£ h )]+Vyo)\S(oxl - 6x 0 (ft) + V£ h )+ 2(2ft> + 3V£„ )] 

3M„liVHrZ-x 2 a ) ' 


Three of these parameters, namely q, s, and r are dependent of the linear control gain y /, . This is 
important toward determination of the flutter boundary, that becomes function of the control 
gain. These expressions substantiate in a strong way the previously mentioned statement, that 
within the present approach, the issues of the control and nature of flutter boundary can be 
determined in an analytic way. 

As a reminder, for steady motion, the equilibrium is stable, in Liapunov’s sense, if the real parts 
of all the roots of the characteristic equation are negative 26 . It is well known that, the analysis of 
the roots of this equation, leads to the Routh-Hurwitz (R-H) conditions, which define the system 
parameter bound for the stability of the considered state of equilibrium. 

The R-H conditions reduce to the inequalities p> 0, q>0, r>0, s>0, and 

pqr-sp 2 -r 2 > 0. For the aeroelastic stability problems, in which the condition 


A, = sp/r+ p 2 j 4 > 0 should be satisfied, the roots of the characteristic equation on the critical 
flutter boundary, SR = 0 , are given by: 

(O l 2 - ±ic , ft), 4 -m± in where i — V— 1 , (24a) 

and 




m = -- 


n> = V_P_ 
4 


n > 0 . 


(24b) 


p 2 r 

The first of equation (24a) reveals that the required condition for the application of Hopf 
bifurcation is fulfilled. 

For sufficient small values of the speed V , all the roots of the characteristic equation are in the 
left hand side half-plane of the complex variable and the zero solution of the system is 
asymptotically stable. The value V = V F for which the two roots of the characteristic equations 
are purely imaginary and the remaining two are complex conjugate and remain also in the left 
hand side half-plane of the complex variable, is critical and corresponds to the critical flutter 
velocity. 

Remembering that the condition of Hopf-Bifurcation is obtained if the characteristic equation has 
a pair of purely imaginary roots, then, in order to identify the “undangerous” portions of the 
stability boundary from the “dangerous” ones is necessary to solve the stability problem for the 
system of equations in state-space form in the critical case of a pair of pure imaginary roots. 


13 



The expression of the boundary of the region of stability 9i for the open/closed loop cross 
sectional wing in a supersonic/hypersonic flow field is defined by the equation: 

* = — nrTr, — -A'O+vfiV- - v ' 2 yk-i)K 2 (vy+2M^c a +2A/^„) 
27MVVfc-^) 

+ Vy[4 + 3xl + 6* n (*„ -l)-6*J 2 -3A^ 2 {3r;[2M„/i« 2 C 0 + (\+yf l XVy + 2M^ k )] 

+ Vy(o[4co + 3(5x1 + 6V£„ - 6x 0 (co + V£„ )f (r a 2 - xl )+ fe [ 2 M^W 2 ^ a + (l + yq ) 

(Vy + 2 )]+ Vy5>\ 4® + 3 ®x 2 a + 6V£„ - 6x 0 (a? + v;„ ®3r „ 2 (Vy + 2M„rf a 
+ 2 M'.timZ" y Ky[4 + 3.v 2 + 6jr 0 (*„ - 1)- 6 Xa ®3 [m ju(l I + « 2 + V', ) 

+ 2C a (Vy + 2 M'.noSt )]+ Vy[SM^ h + 6 M^cox^ h - 3 Af _/tt 0 (V + 4 gJ£, ) 

+ V'(y + 3M^-3M„ W J]}}=0 .(25) 

Notice that this expression is general and includes the relationship between the flutter speed and 
the flutter frequency, parameters evaluated on 9? = 0 in terms of the basic geometrical and flight 
parameters. In the particular case in which the structural damping ratios in plunging and pitching 
are discarded and y => 1 , i.e. no correction for the flight speed range is included, the expressions 

of frequency % F ar| d speed V F on the flutter instability boundary can be obtained as: 


X F ~ 


CO 


Jf 


r a + 0 ~ *0 ) 2 + I ~ 2 Xg (l ~ fo) 

(l + Vfi> a 2 +« [(l-* 0 ) 2 +i] ’ 


bco a 4 x~f V W [(! + Vfi + (1 - JCo X® : Xf - 0 + Vfi ))]- I(i+Vfi) ’ 


(26) 

(27) 


These expressions, constitute the extended counterparts of those obtained by Ashley and 
Zartarian 21 ’ 2:, in the case of absence of the control. As a result, Eqs. (26) and (27) represent the 
dimensionless flutter frequency and flutter speed of the closed-loop system. As a remark, in these 
expressions only the linear control gain is involved. 

In Fig. 3 the influence of the linear control gain i/q on the flutter speed is emphasized. It clearly 
appears that with the increase of the control gain, an increase of the flutter speed is experienced. 
Moreover, toward larger values of t//, the control is more effective and a larger increase in the 
flutter speed is experienced. 

Following the ideas developed by Liapunov 27 - Bautin 16 , the dangerous or undangerous portions 
of the flutter instability boundary can be determined, via determination of the sign of the 
Liapunov’s coefficient. Following Bautin, the system of equations (21) is reduced to the 
canonical form as: 
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(27a) 


^=-c«;+s, 

6fT 

, 

dr 

if=m('-n4:+QMU',4;4:)- 


(27b) 


(27c) 


^=<-m« + e 4 («;,£,£), 

dr 


(27d) 


where m, n and c are determined by Eqs. (24b). In Eqs. (27) Q ( . denote the nonlinear terms 

appearing in the canonical counterpart of Eqs. (21). The expression of the first Liapunov quantity 
is given under a closed form in Ref. 1. For the present case, the first Liapunov quantity is 
expressed in terms of the coefficients as: 

L(V f ) = M! + A m + A \n + A ( 22 ), (28) 

4 c 

where, the terms in the bracket of Eq. (28) are expressed via the coefficients a\ A appearing in 


Eqs. (18), as: 


A kb . ji a 222^2k^2l^2s + ^ j4 a 222&2k^2l&2s )* 


The various components of (28), evaluated on the instability boundary are explicitly expressed 


aO) --LU +d a {4) \y 3 

— A V lx 13 u 222 “ uc |4 a 222 / x 21 ’ 

A n 


(30a) 


a (y) +d flWW 1 

^222 “ A \ U ’ 23 M 222 ~ Uf 24 u 222 7^22 ’ 

A 0 

A U 2 = ~7 (0^73^222 "t" ^“24^222 )^22®l*2 ’ 
An 


(30b) 


(30c) 


aO) (ft +& a^\l 

^122 ” A \ U '13 W 222 ^ oc 14 w 222 H", 

A n 


( Y 2 

) 2 U '22 * 


(30d) 


where the adjoint of a matrix is defined as U=A.kf' and A 0 = j|ot (/ 1| is the determinant of 
the matrix of coefficient a.. . In addition, the parameters a . are: 
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a, 2 = c| 


( 2 ) ( 2 ) 
0 C *d 


a 


o) „(*>i 


3 


al 


a v - 


4" 

<41 

1 Cl 

r 

(3) 

3 

ak 

<*?> 

(1) 

3 

<*, 0) 

«4 W 

.(2) 

3 

,‘ 2) 


a t (,) 

a" 


/V\ 


+ c 


a 4 ] 



• (1) ai 0 


-c 2 a< 2) . 


,P> 

*3 

,( 2 ) 


,(i) 


(31a) 


(31b) 


(31c) 


while the coefficients of the system of Eq. (29) are: 


jb> = 

* 222 


- ik„ + ] 


W 7 ^) 12/l(r a 2 - Z J 

= (*,s« +<s c c„ £ k a -ik« +k] 


(32a) 


«222 =-( 5 S fi + 5 cV' 2 )- 


t + 5. 


Mki+rk 3 (x 0 -i+^) 


V&rl-Xa) A 12/l(r„ 2 -^J 

= — (^5 8 nl + 8 c C nl )+ 5 a A nl {x 0 — 1 + x a ) . (32b) 

Herein the physical and aerodynamic nonlinear parameters as well the nonlinear control terms 
are defined as: 


2 

r a 


8 NL ~ B— 2 / 2 " 1 ’ 

V F Va-Xa) 


M.O+r)*' 

Mpfa-X*)' 


C NL ~ Wz ' 


r a 


(33) 

(34) 

(35) 


VAra-Xa) 

where V F is the flutter speed. Moreover, S A , S s and S c are tracers that, depending on the effect 


that are identified, should be taken as one when the respective effect is taken into account, and 
zero when is discarded. 

Defining fcl-fe u \' and noticing that c>0, the flutter critical boundary is undangerous or 

dangerous, if the following inequalities, L(V f )<0 and L(V' f )>0, are fulfilled, respectively, 
where: 

l(V F )= {d n a^ 2 +^ 14 <3 222 )*21 + (<®23 fl 222 +^24«222 }*22 

+ (<^23«222 + <*24°222 2 + (®|3«222 +^14^222 ^Z^ZZ • (36) 
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In an extended form, the above expression can be cast as: 

Lty f ) = ^C^NL )l(^l $Xa ~ ^14^21 ~~&24 )^22 faliXa ”’**24 )^22^12 ^fi\$Xa ~ ^4)^12^22] 

A^NL {^13 [(*() ~^)Xa "I" r a ] — ^14 t*0 — ^ /ta )}^21 ^23 [(*(> — OZa V a \~ ^24 (*l> “ ^ Xa )}^22 
+ {*2.,[(*0 ~ l )x a + ^]-«24^0 -l + ^a)K2«l2 + {*13 [( *., “ O^a + ^ ]-<* 14 t*o “ 1 + Xa )Fl2«22 } . (37) 

Using the expressions of S NL , A NL and C NL (Eqs. (33) through (35)); from Eq. (37) that is 

■? 

multiplied by V F ~ , paralleling the procedure devised in Refs. 1 and 2, the character of the flutter 
boundary, i.e. undangerous or dangerous, may be expressed, respectively, under the concise 
form as: 

U f 2 <V r 2 , or V f 2 >U r 2 . (38a,b) 

where: 

V r 2 = AjA 2 . (38c) 

In Eq. (38c) the parameter A ] contains the physical and the nonlinear control gain parameter, 
whereas A 2 , includes the aerodynamical nonlinearities. These are expressed as: 


A = £ _ (2 |4 )x\ + (d 23 x a -<* 24 

ta - Xa) 

+ &>Xa -^24^22^2 +(^nXa 

A = y? 2 +7 ^ Pl 3 ho - 1 )Xa+ r a ]~ <*,4 (*0 “ 1 + Xa )}*21 

12 M r a~Xa) 

+ {*23 ho ~ l )Xa + r a I" <*24 (*0 " 1 + Xa )}*22 
+ {*23 [( X 0 ~ l )Xa + r a V <*24 ( X 0 ~ 1 + Xa )j“22«.2 
+ {$ 13 [(jt: 0 — \)Xa r a ] — ^ 14(^0 — 1 + Xa )}*]2<*22 } 


As a very important consequence, in absence of nonlinear control, for B< 0, i.e. for soft 
physical nonlinearities, the expression of V r 2 is negative, and the relation V F > V r 2 is satisfied for 


any flight supersonic Mach number. This implies that in this case, even in the presence of the 
linear control, a subcritical Hopf-Bifurcation is experienced, (see Fig. 9, where the Liapunov 
First Quantity is depicted). On the other hand, for B > 0 , i.e. for hard physical nonlinearities, in 
the presence of the linear control, the transition from the undangerous flutter boundary to the 
dangerous one occurs at higher flight Mach numbers (see, for example Fig. 9). 


Moreover, as mentioned before, if x 0 = 1 , Eq. (37) reduces to the following one: 
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F ) — ^C^NL )l(®-13^a ®I4 )®21 (®23^a ®24 )®22 "*" (®23^o ^24 )®22®I2 (,®\lXa ** 14 )^ 12 ^ 22 ] 

- & A A NL [fe^ - -a u Xa hll + feC - KXa )4 + fee " ^Xa fefe + fe.'a “ KXa fe«22 ] >(40) 


In this case, it clearly appears that a decrease of the influence of the aerodynamic nonlinearities 
on the aeroelastic system is experienced. 


Critical Flutter Boundary: Stability Analysis in 


Presence of Active Control 


As to render clearly the results concerning the critical flutter boundary for open/closed loop 
cross-sectional wings, as displayed in Figs. 6 through 16, some explanations on a couple of 
generic plots. Figs. 4a and 4b, are given next. In Fig. 4a, the intersections of the two curves V F 
and V r , separate the parts of the flutter boundary that are undangerous (for V F <V r ) from the 
dangerous ones, for which the opposite relationship is satisfied, namely V F >V r . Needless to say, 
in such a case a switch of the sign in the first Liapunov quantity occurs, (see Fig. 4b). The graph 
depicting the first Liapunov quantity L(V f ), that defines the portions of the boundary of the 
region of stability, is displayed in Fig. 4b. On the curves that represent the flutter critical 
boundary, the “dangerous” and “undangerous” portions are being indicated for each case by the 
dotted and gross lines, respectively. In the range in which V F < V r , we have the following 
possible situations: i) for V < V F , as time unfolds, a decay of the motion amplitude is 
experienced, which reflects the fact that in this case the subcritical response is involved, ii) an 
ideal condition, for V =V F the center limit cycle occurs, which reflects the fact that in this case a 
periodic orbit is experienced and iii), for the speed parameter V > V F , the response becomes 
supercritical, while for V F > V r the subcritical flutter boundary is experienced. The parameters in 
use for the simulations, unless otherwise specified, are chosen as: 

/j. = \ 00,x a =0.25:07 = 1.2 :r a =0.5;£ a = £„ =0;*,, =0.5;y = 1 ;jc = 1.4;^ =8 S =8 C =1;B = 50. 

The effect of physical nonlinearities on the character of the flutter boundary is carried out in 
terms of the nonlinear parameter B , see Eq. (4). For the condition indicated, the aeroelastic 
system appears to be characterized by subcritical Hopf bifurcation in the upper half plane 
(dangerous flutter boundary) and supercritical in the lower half plane (undangerous flutter 
boundary). The graph depicting the first Liapunov quantity L(V f ) for uncontrolled lifting surface 
and for the cases in which aerodynamic and physical hard nonlinearities are retained, is 
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displayed in Fig. 6. In these simulations both types of nonlinearities, have been considered 
separately and together. It clearly appears that in presence of aerodynamic nonlinearities only, 
the Liapunov quantity is positive for any flight Mach number. This result reflects the fact that 
this type of nonlinearities provides a “catastrophic” character to the flutter boundary, implying 
that a subcritical H-B occurs. On the other hand, in presence of physical hard nonlinearities only 
the opposite situation is experienced. In this sense, at relatively moderate supersonic flight Mach 
numbers, the flutter is benign, while with the increase of the flight Mach number, the flutter 
becomes dangerous. This implies that for larger flight Mach numbers the effects of the 
aerodynamic nonlinearities become prevalent. It is also shown that the neglect of nonlinear 
aerodynamic terms, yields inadvertent results related to the character of the flutter instability 
boundary, especially at high flight Mach numbers, when the aerodynamic nonlinearities become 
detrimental from the point of view of the character of flutter boundary. 

Moreover, as a consequence of these results, if the aerodynamic nonlinearities are discarded 
(5 A = 0) the aeroelastic system features supercriticality or subcriticality for any flight Mach 
number, depending on whether hard ( B > 0) or soft ( B < 0) physical nonlinearities are present, 
respectively. 

The influence of the hard physical and aerodynamic nonlinearities for controlled/uncontrolled 
system is presented in Fig. 7. The dotted and solid lines identify the case in which both physical 
hard ( B = 50 ) and aerodynamical nonlinearities, and aerodynamic nonlinearities only (5 = 0), 
are included, respectively. The control mechanism act in both situations toward the stabilization 
of the system, i.e. to enhance the flutter behavior. Also in this case, the inherent dangerous 
character of the flutter boundary that corresponds to the case when only aerodynamic 
nonlinearities are considered, can be converted into a benign one. 

Figures 6 and 7 show that the soft nonlinearities contribute to render the system unstable and that 
only acting with a nonlinear control the inherent catastrophic type boundary can be converted 
into a benign one (Fig. 6). Moreover, it clearly appears that the linear active control cannot 
change the character of flutter boundary when soft physical nonlinearities are present (Fig. 7). 
From Eq. (28) defining the Liapunov quantity, the supercritical parts of the flutter boundary are 
defined in a closed form. Moreover, by using Eqs. (37), curves associated with V F and V r have 
been plotted in Figs. 10, 12, 14 for 5 = 50 and <5^ = 1 . Each of these graphs present in the plane 
(V, knight) (A = M%, where X = l/(MX a r a ) is the scaling factor), the dangerous and undangerous 
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portions of the flutter boundary for open/closed loop sectional wings. The corresponding 
Liapunov quantity, is depicted in Figs. 11, 13, 15, in the plane (L, A-night)- Moreover, in Figs. 10 
through 16 aerodynamic and physical hard nonlinearities have been included <5^=1, <5 S = 1 ; 

B = 50 . With dotted lines are marked the values of the flight Mach number where the transition 
from the benign to the catastrophic character of the flutter boundary occurs. As a consequence, a 
complete picture of the variation of dangerous and benign parts of the flutter boundary has been 
provided. In the plots depicting the Liapunov quantity, the gain parameters y/ : ( i=l,2 ), help to 

understand the potentiality of the linear and nonlinear active control to enhance the flutter 
instability behavior and convert the dangerous flutter boundary into an undangerous one. 

Notice that, the linear active control is able to act also toward the enhancement of the flutter 
boundary. This can be readily seen either from Eq. (27) and from Figs. 3 and 5 where for if/ ] * 0 , 
V F was depicted, and the beneficial influence on the increase of the flutter speed revealed. 

The dangerous and undangerous portions of the flutter boundary for open/closed loop sectional 
wings via the use of the two different aerodynamic theories (PTA and SWT) are presented in 
Fig. 16. In addition to what was stated before about these two theories, from this plot and from 
Table 1, it can be revealed that the PTA gives more conservative results (less then 3% difference 
at the transition flight Mach number), as compared with the SWT, in the sense that within PTA 
the transition from the benign to catastrophic flutter boundary appears at slightly lower flight 
Mach numbers as compared with those featured by SWT. This conclusion was expected by 
comparing the cubic nonlinearities in the two pressure expressions (see Eqs. (10) and (12)), 
wherefrom it clearly appears that the one within the PTA is larger than that in SWT. 

Conclusions 

It was shown that in some circumstance the hard physical and aerodynamic nonlinearities 
contribute in different manners to the determination of the “undangerous” or “dangerous” 
character of the flutter critical boundary, in the sense that the “hard” nonlinearities render the 
flutter boundary benign, whereas the “soft” ones contribute to its catastrophic character. In 
addition, at high flight Mach numbers the aerodynamic nonlinearities contribute invariably to the 
“dangerous” character of the flutter boundary. This implies that, with the increase of the 
hypersonic flight speed, when the aerodynamic nonlinearities become prevalent, the flutter 
boundary will become catastrophic, irrespective of the presence of hard physical nonlinearities. 
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On the other hand, soft physical nonlinearities (B < 0) contribute in the same sense, as the 
aerodynamic nonlinearities, to the dangerous character of the flutter boundary. It was also shown 
that, in the case when due to the nature of the involved nonlinearities, when the catastrophic 
aeroelastic failure can occur, an active control can be used as to convert the flutter boundary in a 
benign one or to shift the transition between these two states toward larger flight Mach numbers. 
Numerical simulations illustrate that: i) the increase of the values of the active control gains 
yields an increase of the “undangerous” portions of the flutter instability boundary, ii) with the 
increase of the supersonic flight Mach numbers that results in an increase of aerodynamic 
nonlinearities, a decrease of the “undangerous” portions of the flutter instability boundary is 
experienced, iii) in the cases in which catastrophic aeroelastic failures can be experienced, the 
active control capabilities can render the flutter boundary safe, in the sense that its crossing does 
not yields a catastrophic failure, and iv) a comparisons of the results provided when PTA and the 
SWT are used, reveal that the PTA gives more conservative results. 

As clearly appears from this paper, the issue of generating the active control moment was not 
addressed. It is in the authors’ believe that this can be produced via a device operating similarly 
to a spring, whose linear and nonlinear characteristics can be controlled. 

To the best of the authors’ knowledge, this is the first paper where the problem of the character 
of the flutter boundary of cross-sectional wing was approached in a so general context, in the 
sense that, the implications of the various nonlinearities considered in conjunction with the 
presence of a control capabilities have been considered together as to analyze this problem, and 
determine also the combination of those nonlinearities and parameters that can yield a 
conversion of the catastrophic type of flutter into a benign one. Needless to say, the concept and 
methodology presented here can be extended as to treat more complex associated problems 
involving 3-D straight and swept aircraft wings. 
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Appendix A 

The coefficients of the aeroelastic governing system represented in state-space form, Eqs. (21): 
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Figure Captions 

Figs, la and lb: Character of the flutter boundary in the terms of LCOs amplitudes - Vp flutter 
speed of the linear system; H-B: Hopf-Bifurcation. 

Fig. 2: Comparison of the predictions of the flutter speed vs. the flight Mach number when using 
Piston Theory Aerodynamics (PTA), of aerodynamics based on Shock Waves Theory (SWT), 
and of exact unsteady supersonic aerodynamics, 

(fi =lQO;x a =0.25;<y =1.2;r 0 =0.5;£ a =C h = 0;* 0 =0.5). 

Fig. 3: Flutter speed vs. flight Mach number. Influence of the linear control gain t //, . 

Figs. 4a and 4.b: Generic representation of the dangerous and undangerous portion of the flutter 
critical boundary. The upper plot (4. a) is in the (V - M^ght) plane, while the bottom one (4.b) is 
in the (L - Mnight) plane; L < 0 correspond to V r > Vp and vice versa. 

Fig. 5: Evaluation of the flutter speed for selected flight Mach numbers vs. linear control gain 
W 

Fig. 6: Influence of the physical and aerodynamical nonlinearities on the First Liapunov 
Quantity L for the uncontrolled cross-sectional wing. 

Fig. 7: Influence of the physical nonlinearities on the First Liapunov Quantity L in the 
presence/absence of linear and nonlinear active control gains. Aerodynamic nonlinearities 
retained. Flutter character change for system encompassing physical hard nonlinearity. 

Fig. 8: Conversion from subcritical to supercritical H-B via nonlinear active control, for system 
encompassing physical soft nonlinearity ( B = — 10). 

Fig. 9: Subcritical H-B for system encompassing physical soft nonlinearity (5 = -10), linear 


active control included. 



Fig. 10: Undangerous and dangerous portions of the flutter critical boundary in the presence of 
linear control. 

Fig. 11: Influence of the linear control on the First Liapunov Quantity L. 

Fig. 12: Undangerous and dangerous portions of the flutter critical boundary in the presence of 
nonlinear control. 

Fig. 13: Influence of the nonlinear control on the First Liapunov Quantity L. 

Fig. 14: Undangerous and dangerous portions of the flutter critical boundary in the presence of 
linear and nonlinear controls. 

Fig. 15: Influence of the linear and nonlinear controls on the First Liapunov Quantity L. 

Fig. 16: Undangerous and dangerous portions of the flutter critical boundary with and without 
control in presence of both nonlinearities. Comparison between SWT and PTA. 
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